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Japanese laser interferometric gravitational wave detectors, TAMA300 and LISM, performed a 
coincident observation during 2001. We perform a coincidence analysis to search for inspiraling 
compact binaries. The length of data used for the coincidence analysis is 275 hours when both 
TAMA300 and LISM detectors are operated simultaneously. TAMA300 and LISM data are analyzed 
by matched filtering, and candidates for gravitational wave events are obtained. If there is a true 
gravitational wave signal, it should appear in both data of detectors with consistent waveforms 
characterized by masses of stars, amplitude of the signal, the coalescence time and so on. We 
introduce a set of coincidence conditions of the parameters, and search for coincident events. This 
procedure reduces the number of fake events considerably, by a factor ~ 10 -4 compared with the 
number of fake events in single detector analysis. We find that the number of events after imposing 
the coincidence conditions is consistent with the number of accidental coincidences produced purely 
by noise. We thus find no evidence of gravitational wave signals. We obtain an upper limit of 0.046 
/hours (CL = 90%) to the Galactic event rate within lkpc from the Earth. The method used in 
this paper can be applied straightforwardly to the case of coincidence observations with more than 
two detectors with arbitrary arm directions. 
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I. INTRODUCTION 

In the past several years, there has been substantial progress in gravitational wave detection experi- 
ments by the ground-based laser interferometers, LIGO[l|, VIRGO0, GEO6OO0, and TAMA300[1@. 
The observation of gravitational waves will not only be a powerful tool to test general relativity, but also 
be a new tool to investigate various unsolved astronomical problems and to find new objects which were 
not seen by other observational methods. 

The Japanese two laser interferometers, TAMA300 and LISM, performed a coincident observation 
during August 1 and September 20, 2001 (JST). Both detectors showed sufficient stability that was 
acceptable for an analysis to search for gravitational wave signals. Given the sufficient amount of data, 
it was a very good opportunity to perform a coincidence analysis with real interferometers' data. 

There were several works to search for gravitational waves using interferometeric data. A coincidence 
analysis searching for generic gravitational wave bursts in a pair of laser interferometers has been reported 
in |g| . Allen et al.0 analyzed LIGO 40m data and obtained an upper limit of 0.5/hour (CL = 90%) on 
the Galactic event rate of the coalescence of neutron star binaries with mass between 1M Q and 3M©. 
Tagoshi et al.Q analyzed TAMA300 data taken during 1999 and obtained an upper limit of 0.59/hour 
(CL = 90%) on the event rate of inspirals of compact binaries with mass between 0.3M Q and 10M© and 
with signal-to-noise ratio greater than 7.2. Very recently, an analysis using the first scientific data of the 
three LIGO detectors was reported 0, and an upper limit of 1.7 x 10 2 per year per Milky Way Equivalent 
Galaxy is reported. Recently, International Gravitational Event Collaboration (IGEC) of bar detectors 
reported their analysis using four years of data to search for gravitational wave bursts [13. They found 
that the event rate they obtained was consistent with the background of the detectors' noise. 

In the matched filtering analysis using real data of single laser interferometer (e.g. 0, 0), many fake 
events were produced by non-Gaussian and non-stationary noise. In order to remove such fake events, 
it is useful to perform coincidence analysis between two or more independent detectors. Furthermore, 
coincidence analysis is indispensable to confirm the detection of gravitational waves when candidates for 
real gravitational wave signals are obtained. The purpose of this paper is to perform coincidence analysis 
using the real data of TAMA300 and LISM. 

We consider gravitational waves from inspiraling compact binaries, comprized of neutron stars or black 
holes. They are consider to be one of the most promising sources for ground based laser interferometers. 
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Since the waveforms of the inspiraling compact binaries are known accurately, we employ the matched 
filtering by using the theoretical waveforms as templates. Matched filtering is the optimal detection 
strategy in the case of stationary and Gaussian noise of detector. However, since the detectors' noise is 
not stationary and Gaussian in the real laser interferomters, we introduce x 2 selection method to the 
matched filtering. 

We analyze the data from each detector by matched filtering which produces event lists. Each event 
is characterized by the time of coalescence, masses of the two stars, and the amplitude of the signal. If 
there is a real gravitational wave event, there must be an event in each of the event list with consistent 
values of parameters. We define a set of coincidence conditions to search for coincident events in the two 
detectors. We find that we can reduce the number of events to about 10 -4 times the original number. 
The coincidence conditions are tested by injecting the simulated inspiraling waves into the data and by 
checking the detection efficiency. We find that the detection efficiency is not affected significantly by 
imposing the coincidence conditions. 

We estimate the number of coincident events produced accidentally by the instrumental noise. By using 
a technique of shifting the time series of data artificially, we find that the number of events survived after 
imposing the coincidence conditions is consistent with the number of accidental coincidences produced 
purely by noise. 

We propose a method to set an upper limit to the real event rate using results of the coincidence 
analysis. In the case of TAMA300 and LISM, we obtain an upper limit of the event rate as 0.046/hour 
(CL = 90%) for inspiraling compact binaries with mass between IMq and 2M Q which are located within 
lkpc from the Earth. In this case, since TAMA300 is much more sensitive than LISM, the upper limit 
obtained from the coincidence analysis is less stringent than that obtained from the TAMA300 single 
detector data analysis. This is because the detection efficiency in the coincidence analysis is determined 
by the sensitivity of LISM. Thus, the upper limit obtained here is not the optimal one which we could 
obtain using the TAMA300 data taken during 2001. 

The method to set an upper limit to the event rate proposed here can be extended straightforwardly 
to the case of a coincidence analysis for a network of interferometric gravitational wave detectors. 

This paper is organized as follows. In Section[nJ we briefly describe the TAMA300 and LISM detectors. 
In Section ITTT1 we discuss a method of matched filtering search used for TAMA300 and LISM data. In 
Section IIVI the results of the matched filtering search for each detector are shown. In section [V] we 
discuss a method of the coincidence analysis using the results of single-detector searches, and the result 
of the coincidence analysis is shown. We also derive the upper limit to the event rate in Section IVII 
Section [VIII is devoted to summary. In Appendix A, we discuss a % 2 veto method to distinguish between 
real events and fake events produced by non-Gaussian noise. In Appendix B, we examine a different 
choice of At ( the length of duration to find local maximum of matched filtering output ) for comparison. 
In Appendix C, we discuss a sidereal time distribution of coincidence events. In Appendix D, we review 
a method to estimate the errors in the parameters due to noise using the Fisher matrix. 

Throughout this paper, the Fourier transform of a function h(t) is denoted by h(f), which is defined 

by 

poo 

h(f) = / dt e^hit). (1) 



II. DETECTOR 
A. TAMA300 

TAMA300 is a Fabry-Perot-Michelson interferometer with the baseline length of 300m located at the 
National Astronomical Observatory of Japan in Mitaka, Tokyo (35.68°N, 139.54°E) (See Table HJ). The 
detector's arm orientation (the direction of the bisector of two arms) measured counterclockwise from 
East is 225°. The details of TAMA300 detector configuration can be found in @. The TAMA300 detector 
became ready to operate in the summer 1999 Most of the designed system (except power recycling) 
were installed by the that time. First data taking was performed as a test during August 1999 (DTI). In 
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frequency [Hz] 

FIG. 1: The strain equivalent noise spectra of TAMA300 and LISM on September 3, 2001. 

September 1999, three days observation (DT2) was carried out, and the first search for gravitational waves 
from inspiraling compact binaries was performed [8|. Since then, TAMA300 has been performing several 
observations. In August 2000, an observation (DT4) was performed for two weeks and 160 hours of data 
were taken which are described in detail in Q. From March 2nd to March 8th, 2001, TAMA300 performed 
an observation (DT5) and 111 hours of data were taken. After improvements of the sensitivity, TAMA300 
had carried out a long observation (DT6) from August 1st to September 20th, 2001. The length of data 
taken was about 1100 hours. The best strain equivalent sensitivity was about fe~5x 10~ 21 /\/Hz around 
800Hz at DT6. From August 31th to September 2nd, 2002, TAMA300 performed a short observation 
(DT7) and 24 hours of data were taken. From February 14th and April 15th 2003, TAMA300 performed 
an observation (DT8) for two months, and 1158 hours of data were taken. Most recently, from November 
28th 2003 to 10th January, 2004, TAMA300 performed an observation (DT9) and 557 hours of data were 
taken. The observation history of TAMA300 is summarized in Table ITT1 

In this paper, we use the DT6 data taken from September 2nd to 17th, 2001 when LISM was also in 
good condition. The amount of data available for the coincidence analysis is 275 hours in total. Typical 
one-sided noise power spectra of TAMA300 and LISM during this observation are shown in Fig. ^ 



B. LISM 

LISM is a laser interferometer gravitational wave antenna with arm length of 20m, located in the 
Kamioka mine (36.25°N, 137.18°E), 219.02km west of Tokyo. The detector's arm orientation is 165° 
measured counter clockwise from East. The LISM antenna was originally developed as a prototype 
detector from 1991 to 1998 at the National Astronomical Observatory of Japan, in Mitaka, Tokyo, to 
demonstrate advanced technologies [TT] . In 1999, it was moved to the Kamioka mine in order to perform 
long-term, stable observations. Details of the LISM detector is found in 

The laboratory site is 1000m underground in the Kamioka mine. The primary benefit of this location 
is extremely low seismic noise level except artificial seismic excitations. Furthermore, much smaller 
environmental variations at this underground site are beneficial to stable operation of a high-sensitivity 
laser interferometer. The optical configuration is the Locked Fabry-Perot interferometer. The finesse of 
each arm cavity was about 25000 to have a cavity pole frequency of 150Hz. The main interferometer was 
illuminated by a Nd:YAG laser yielding 700mW of output power, and the detector sensitivity spectrum 
was shot-noise limited at frequencies above about 1kHz. 
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TAMA300 (DT6) 




LISM 


Interferometer type 


Fabry- Perot-Michelson 




Locked Fabry-Perot 


Base length 


300m 




20m 


Finesse of main cavity 


500 




25000 


Laser Source 


NdYAG, 10W 




Nd:YAG, 700mW 


Best sensitivity in strain h [1/VHz] 


5 x 10 -21 




6.5 x 10~ 20 


Location and arm orientation 


35.68°N, 139.54°E, 225° 


36.25°N, 137.18°E, 165° 


Maximum delay of signal arrival time 


0.73msec 




Operation period 


Aug.l - Sept.20, 2001 


Aug 


.1 - 23, Sept.3 - 17, 2001 


Observation time 


1038 hours 




786 hours 


Operation rate 


87% 




91% 


Simultaneous observation 


709 hours 




Data used for coincidence analysis 


275 hours 





TABLE I: Summary of the observation in August and September 2001 by TAMA300 and LISM 





Year 


period 


obsevation time [hours] 


Topics 


DTI 


1999 


6-7 Aug. 


11 


Total detector system check 










and Calibration test 


DT2 


1999 


17-20 Sept. 


31 


First event search 


DT3 


2000 


20-23 April 


13 


Sensitivity improved 


DT4 


2000 


21 Aug.-4 Sept. 


167 


100 hours observation 


DT5 


2001 


2-10 Mar. 


111 


Full time observation 


DT6 


2001 


1 Aug.-20 Sept. 


1038 


1000 hours observation 










and coincident observation with LISM 


DT7 


2002 


31 Aug.-2 Sept. 


25 


Power recycling installed (Full configuration) 


DT8 


2003 


14 Feb.- 14 April 


1158 


Coincident observation with LIGO 



DT9 2003 - 2004 28 Nov.- 10 Jan. 557 Full automatic operation 

and Partial coincident observation 
with LIGO and GEO600 



TABLE II: Observation history of TAMA300 



The operation of LISM was started in early 2000, and has repeatedly been tested and improved since. 
The data used in this analysis were taken in the observations between August 1st and 23th and between 
September 3rd and 17th, 2001. The total length of data is 780 hours. The first half of the period was 
in a test-run and some improvements were made after that. The data from the second half were of good 
quality to be suitable for a gravitational wave event search, so 323 hours of data for the latter half was 
dedicated for this analysis. The best sensitivity during this period was about h <~ 6.5 x 10 _20 /\/Hz 
around 800Hz. 
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III. ANALYSIS METHOD 
A. Matched filtering 

To search for gravitational waves emitted from inspiraling compact binaries, we use the matched 
filtering. In this method, cross-correlation between observed data and predicted waveforms are calculated 
to find signals and to estimate binary's parameters. When the noise of a detector is Gaussian and 
stationary, the matched filtering is the optimal detection strategy in the sense that it gives the maximum 
detection probability for a given false alarm probability. 

We use restricted post-Newtonian waveforms as templates: the phase evolution is calculated to 2.5 
post-Newtonian order, and the amplitude evolution is calculated to the Newtonian quadrupole order. 
The effects of spin angular momentum are not taken into account here. The filters are constructed in 
Fourier domain by the stationary phase approximation p| of the post-Newtonian waveforms 0]. We 
introduce the normalized templates h c and h s which are given in the frequency domain for / > by 

h c = Nf- 7 / 6 exp(iV(f)), (2) 
h s = iNf- r / e exp(iy(f)), (3) 



where 



*(/) = 2nft c - l + JL-^GMfc- 3 )- 5 / 3 



1 + K^TT + (irGMfc- 3 ) 2 / 3 - 16n(nM fc~ 3 ) 
9 V 84 J 



( 15293365 27145 
V 508032 + 504 



3085 2 \ , 3x4/3 tt / 38645 \. „ rJ , 3, 5/3 



(4) 



where / is the frequency of gravitational waves, t c is the coalescence time, M = mi +m2, i] = mxmijM 2 ! , 
and mi and ttt-2 are the masses of binary stars. For / < 0, they are given by h c / s (f) = h*, s (—f), where 
the asterisk denotes the complex conjugation. The normalization factor N is defined such that h c and 
h s satisfy 

(h c ,h c ) = l, (h s ,h s ) = l, (5) 

where 

M - 2 l-J^Jm- (6) 

S n (f) is the strain equivalent one-sided noise power spectrum density of a detector. We note that, for h c 
and h s calculated by the stationary phase approximation, we have (h c , h s ) = 0. 
In the matched filtering, we define the filtered output by 

p = (s, h c cos(4> c ) + h s sm((p c ))) (7) 

where s(t) is the signal from a detector and <p c is the phase of the template waveform. For a given interval 
of t c , we maximize p over the parameters t c ,M, r\ and <f> c . The filtered output maximized over (j) c is given 

by 



p = y/( a ,h c )* + ( S ,h s )* = p. (8) 

The square of the filtered output, p 2 , has an expectation value 2 in the presence of only Gaussian noise 
in the data s(t). Thus, we define the signal-to-noise ratio, SNR, by p/%/2. 

Matched filtering is the optimal detection strategy in the case of stationary and Gaussian noise of 
detector. However, since the detectors' noise is not stationary and Gaussian in the case of real laser 
interferomters, we introduce x 2 method to the matched filtering in order to discriminate such noise from 
real gravitational wave signals. We describe details of x 2 method in Appendix lAl 
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B. Algorithm of the matched filtering analysis 



In this subsection, we describe a method to analyze time sequential data from the detectors by matched 
filtering. 

First, we introduce, "a continuously locked segment". The TAMA300 and LISM observations were 
sometimes interrupted by the failure of the detectors to function normally, which are usually called 
"unlock" of the detectors, or were interrupted manually in order to make adjustments to the instruments. 
A continuously locked segment is a period in which the detector is continuously operated without any 
interruptions and the data is taken with no dead time. In the analysis of this paper, we treat only the 
data in such locked segments. 

The time sequential voltage data of a continuously locked segment are divided into small subsets of 
data with length of 52.4288 seconds (= sampling interval [s] x number of samples = (5 x 10~ 5 ) x 2 20 [s]). 
Each subset of data has overlapping portions with adjacent subsets for 4.0 seconds in order not to lose 
signals which lie across borders of two adjacent subsets. The data of a subset are Fourier transformed 
into frequency domain and are multiplied by the transfer function to transform into strain equivalent 
data. The resulting subset of data is the signal of the detector in the frequency domain, s(/), used in 
the matched filtering. 

The power spectrum density of noise S n (f) is basically evaluated in a subset of data neighbouring to 
each s(t) except for the cases below. On estimating the noise power spectrum, S n (f), we do not use the 
data contaminated by transient burst noise. For this purpose, we evaluate the fluctuation of the noise 
power defined by 



for each set of data with length of 65.6 seconds which composes one file of stored data. We also calculate 
the average of p, (p), within each continuously locked segment. For each s(f), we then apply the following 
criterion. If a subset of data in the neighborhood of s(/) lies entirely in one of the files, we examine the 
value of p of the file, and if it deviates from the average (p) for more than 2dB, i.e., p > 1.26(p), we do 
not use that subset of data for evaluating the power spectrum and move to the neighboring subset. If a 
neighboring subset lies over two files, we examine the values of p of the two files, and if either of them 
exceeds the 2dB level, we use neither of them. If a neiboring subset such that a file (or two consecutive 
files) that contains it has p < 1.26(p) is found, the subset is divided into 8 pieces and the S n (f) is 
evaluated by taking the average of them. If the fluctuations of p are too large, and we cannot find files 
with the values of p within 2dB of the average within the locked segment, we use the power spectrum 
which is evaluated by taking the average of all the data in the corresponding locked segment. 

In order to take the maximization of p in Eq.JSJ over the mass parameters, we introduce a grid in the 
mass parameter space. Each grid poin t defines the mass parameters which characterize a template. We 
adopt the algorithm introduced in [15J to define the grid point in the mass parameter space. The distance 
between the grid points is determined so as not to lose more than 3 % of signal-to-noise ratio due to 
mismatch between actual mass parameters and those at grid points. Accordingly, the mass parameter 
space depends on the power spectrum of noise. In order to take into account of the changes in the noise 
power spectrum with time, we use different mass parameter spaces for different locked segments. For 
each locked segment, the averaged power spectrum of noise is used to determine the grid spacing in the 
mass parameter space. 

We consider the mass of each component star in the range IMq < mi, m-i < 2Mq. This mass range 
is chosen so that it covers the most probable mass of a neutron star, ~ 1.4M©. 

With s(f), S n (f) and a template on each grid point of the mass parameter space, we calculate p in Eq. 
©. For each interval At = 25.6 msec, we search for t c at which the local maximum of p is realized. If 
the p thus obtained is greater than a pre-determined value p m , we calculate the value of \ 2 as discussed 
in Appendix^ We adopt p m = 7 in this paper. Choosing a too large p m results in missing actual events 
from the data, while a too small p m requires too much computational time. The same computation is 
done for all the mass parameters on each grid point. 

Finally, for each interval of the coalescence time with length At = 25.6 msec, we search for t c , M, rj 




-1/2 



(9) 
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FIG. 2: Scatter plots (p, VX 2 ) of the events of TAMA300. 



which realize the local maximum of p. Each maximum is considered a event. The value of t c , /?, % 2 , M, 
r/ of each event are recorded in event lists. 



In this section, we show the result of the independent analysis for each detector. 

Our analysis is carried out with 9 Alpha computers and also with 12 Pentium4 computers at Osaka 
University. The matched filtering codes are paralleled by the MPI library. Among the data from Septem- 
ber 3rd to 17th, 2001, TAMA300 has 292.4 hours of data after removing unlocked periods. We also 
removed the data segments of lengths less than 10 minutes. The total length of data is 287.6 hours. 
LISM has 323.0 hours of data after removing unlocked periods. After removing the data segments less 
than 10 minutes, the total length of data is 322.6 hours. 

The scatter plots of (p, \J~x*) of the events are shown in Figs. and |3| We discriminate the non- 
Gaussian noise from real gravitational wave signals by setting the threshold to the value of /VVa ( see 
Appendix lA"|l . In Figs. 21and[S] > we show the number of events for bins of pj Va^ 

Although the main topic of this paper is to perform a coincidence analysis, for the purpose of comparison 
between a single-detector analysis and a coincidence analysis, we evaluate the upper limit to the event 
rate which is derived from an analysis independently done for each detector. The upper limit to the 
Galactic event rate is calculated by 



where N is the upper limit to the average number of events with pj greater than a pre-determined 
threshold, T is the total length of data [hours] and e is the detection probability. 

To examine the detection probability of the Galactic neutron star binary events, we use a model of the 
distribution of neutron star binaries in our Galaxy which is given by [l6| 



IV. RESULTS OF MATCHED FILTERING SEARCH 




(10) 




dN = e~ R2/2R °e- z/h *RdRdZ, 



(11) 



where R is Galactic radius, Rq = 4.8 kpc, Z is height off the Galactic plane and h z = 1 kpc is the scale 
height. We assume that the mass distribution is uniform between 1M Q and 2Mq. We also assume uniform 
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FIG. 3: The same figure as FigEJbut for LISM. 



TAMA300 

7 1 1 1 1 1 1 



6 



c 




p/(x 2 ) 1 ' 2 

FIG. 4: Histogram of the number of events of TAMA300 in terms of pj 1 \fx* ■ 

distributions for the inclination angle and the phase of an event. With these distribution functions, we 
perform a Monte Carlo simulation. The simulated gravitational wave events are injected into the data of 
each detector for about every 15 minutes. We perform a search using the same code used in our matched 
filter analysis, and evaluate the detection probability for each pj \f~x* threshold. The result for TAMA300 
is shown in Fig. [S] 

For the case of LISM, since LISM's sensitivity is not good enough to observe events in all of the Galaxy, 
we only evaluate the detection probability of nearby events within lkpc. The result is shown in Fig. [7| 

The threshold of p/yx* for each of the analysis is determined by the fake event rate. We set the 
fake event rate to be 2.0 [1/yr]. We approximate the distribution of pj ' \px* in each of Figs. 01 and by 
an exponential function and extrapolate it to large p/va ■ We assume that this function describes the 
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p/(x 2 ) 1/2 

FIG. 5: The same figure as FigEJbut for LISM. 



Galactic event efficiency 

0.7 1 1 1 




5 10 15 20 25 

p/(X 2 ) 1 ' 2 threshold 

FIG. 6: Galactic event detection efficiency of TAMA300. The error bars shows the la error of the simulation. 



background fake event distribution. 

For the TAMA300 case, the fake event rate N hg /T = 2.0 [l/yr]= 0.00023 [1/hour] gives the total 
number of expected fake events as A^g = 0.066. This determines the threshold to be p/ ' yfx* — 14.8. 
With this threshold, we obtain the detection probability, e = 0.263, from Fig. On the other hand, the 
number of observed events with pj \/x* greater than the threshold is iV bs = 0. Using Bayesian statistics, 
and assuming uniform prior probability for the real event rate and the Poisson distributions for real and 
background events, we estimate the expected number of real events A" which pj \fy^ is greater than the 
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FIG. 7: The detection efficiency of TAMA300 and LISM for nearby events within lkpc. The error bars show the 
lcr error of the simulation. 



threshold with a given confidence level (CL). Namely, it can be evaluated from the equation [l7| . 



(N+N bg ) Tpn=N obs (N+N bg ) n 







-JV h 



En— jv 
n=0 



1 - CL. 



(12) 



Using this formula, we obtain the upper limit to the expected number of real events to be 2.30 with 
90% CL. Then, using the length of data T = 287.6 hours, we obtain the upper limit of the event rate as 
R 90% = 0.030 [1/hour] (CL = 90%). 

For the LISM detector, we only evaluate the upper limit to nearby events within lkpc. We set the 



threshold pj \J x 2 = 14.6, corresponding to the number of expected fake events Nb s — 0.074 which realizes 



the fake event rate N^ g /T = 2.0[l/yr]. The number of observed events with p/vx 2 greater than the 
threshold is iV bs = 0. Thus, the upper limit to the expected number of real events is again 2.30 with 90% 
CL. The detection probability is given from Fig[7]as e = 0.042. The length of data is T — 322.6 hours. 
Using these numbers, we obtain the upper limit to the nearby event rate as 0.17 [1/hour] with 90% CL. 



The results of matched filtering analysis for TAMA300 and LISM are summarized in Table IITT 



threshold 


N 


Detection efficiency Length of data 


Upper limit (90% CL) 


TAMA300 14.8 
LISM 14.6 


2.30 (90% CL) 
2.30 (90% CL) 


0.263 287.6 [hours] 
0.042 322.6 [hours] 0.17 


0.030 [1/hour] 
[1/hour] (for nearby events) 



TABLE III: Results of matched filtering analysis for TAMA300 and LISM 



12 



V. COINCIDENCE ANALYSIS 



A. Method 



In the previous section, we obtained event lists for TAMA300 and LISM. Each event is characterized by 
t c , -A4, r], p, and x 2 > where M is the chirp mass (= Mr/ 3 ' 5 ). True gravitational wave events will appear in 
both event lists with different values of these parameters according to the detectors' noise, the difference 
in the detectors' locations and thier arm orientations, and the discreteness of the template space. In this 
section, we evaluate the difference of the parameters real events have. 



Time selection: 

The distance between the TAMA300 site and the LISM site is 219.02km. Therefore, the maximum 
delay of the arrival time of gravitational wave signals is Aidist = 0.73057 msec. The allowed difference in 
t c is set as follows. If the parameter, t CJAMA and t CtL1SM , of an event satisfy 

|^C,TAMA tc, LISM | ^ Aidist At no [ SC: (13) 

the event is recorded in the list as a candidate for real events. We estimate errors in t c due to noise Ai no ; se 
by using the Fisher information matrix (see Appendix iDl for a detailed discussion). We denote the 1 a 
value of the error of t c by At Ct i for i = TAMA or LISM. We determine Ai no i se as Ai no ; sc = a w x At c 

where At c = y/ Ai 2 iTAMA + Ai 2 LISM . The parameter a w is to be determined in such a way that it is small 

enough to exclude accidental coincidence events effectively but is large enough to make the probability 
for missing a real event sufficiently small. 

In this paper, we adopt a w = 3.29 which corresponds to 0.1% probability of losing real signals if the 
noise are Gaussian and if both detectors are located at the same site. Although it may be possible to tune 
the value of o~ w to obtain a better detection efficiency while keeping the fake event rate low enough, we do 
not bother to do so. Instead, we check whether we have a reasonable detection efficiency by this choice. 
To check the detection efficiency is important in any case, since the At c determined above assumes a large 
signal amplitude in the presence of Gaussian noise. The actual detection efficiency might be different 
from what we expected. 



Mass selection: 

In the same way as for t Cl errors in the values of M. and r\ due to detector noise, AA^ no isc and A?7 no i se , 
are estimated by using the Fisher matrix. We denote the 1 a values of errors in M. and r\ by AM.i and Ar/i, 
respectively. We set AM aoise = a w ^J {AM)% AMA + {AM)l lSM and Ar/ noiso = a wy /{An)l AMA + {An)l lSM , 
and adopt a w = 3.29 as in the case of t c . 

When the amplitude of a signal is very large, errors due to detector noise become small since they 
are inversely proportional to p, and errors due to the discreteness of the mass parameter space become 
dominant. We denote the latter errors by AA'lmcsh and Ary mcs h. They are determined from the maximum 
difference in the neighbouring mesh points in the mass parameter space. 

By taking account of the above two effects, we choose the allowable difference in the mass parameters 

as 

|.Mtama - AIlismI < AA4 noise + AM mesh : 

(14) 

Ittama - ??lism| < A77 no i SO + A?7 mcs h. (15) 



Amplitude selection: 

Since the two detectors have different sensitivities, signal-to-noise ratios of an observed gravitational 
wave signal will be different for the two detectors. Further, since their arm orientations are different, the 
signal-to-noise ratios will differ even if they have the same noise power spectrum. 

We express the allowable difference in p TAM A and p LISM as 

^sens ^arm ^noisc — k}§ ( j — ^sens ~t~ ^arm ~t~ ^noisc- (^^) 

^ Plism ' 
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FIG. 8: Relative detection efficiency of the coincidence analysis compared to the single-detector efficiency of LISM 
as a function of the parameter a m used for the coincidence criterion. The dot-dashed line is the efficiency after 
the time selection, the dashed line is the efficiency after the time and mass selection, and the solid line is the 
efficiency after the time-mass-amplitude selection. 



Here, (5 SO ns is due to the difference in S n , 

<5 scns ^ log [f / ^ df) 1/2 /( f f ~ V3 df) 1/2 ] (17) 
J S n TAMA Sn LISM (/) ' J' 

and 5 aTm is due to the difference in the arm orientations, and 5 no i SC is due to detector noise. The value of 
<5noisc is evaluated by the Fisher matrix in the same way as t c and masses. 

The value of S sens is determined for each event individually from the noise power spectrum used in the 
matched filtering. 5 arm is evaluated by a Monte Carlo simulation as follows. We assume that the two 
detectors have the same noise power spectrum, and generate the waveforms of Galactic events randomly. 
We then evaluate p of all the events detected by each detector, and determine the value of <5 arm in such 
a way that for more than 99.9 % of events, we have | log(/9 TAMA //9 LISM )| < S aTm . This gives S aim = 1.60. 



B. Detection efficiency and the parameter windows 



Here, we discuss the detection efficiency of our coincidence analysis. In particular, we examine the 
validity of the choice a w = 3.29 made in the previous section. 

For the Galactic event simulation discussed in Section IIVI the detection efficiencies of TAMA300 and 
LISM for the threshold pj \px^ > 7 are 99% and 24%, respectively. The detection efficiency of the 
coincidence analysis is dominated by the LISM's efficiency. Thus we define the detection efficiency for 
the coincidence analysis, as the fraction of LISM events which fulfill the coincidence criteria. The result 
is shown in Fig. [8] We find that more than 94 % of LISM events can be detected if we set a w > 3. Thus 
with a w = 3.29, we have a reasonably high detection efficiency. 

If we adopt a larger value of a w , we obtain a higher detection efficiency, but the number of fake events 
will also increase, and vise versa for a smaller value of a w . Then, one may tune the value of a w so that 
it gives the most stringent upper limit to the event rate. However, since we cannot expect any drastic 
improvement by such an optimization, we adopt a w = 3.29 in this paper for the sake of simplicity of the 
analysis. 
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FIG. 9: (ptama ,Plism) scatter plots. The crosses (+) are the events survived after the time selection, and the 
circled crosses (©) are the events survived after the time, mass and amplitude selections. 



C. Results 



In this subsection we discuss the results of the coincidence analysis. The length of data used for the 
coincidence analysis is 275.3 hours when both TAMA300 and LISM detector were operated simultane- 
ously. 

As a result of independent matched filtering searches, we obtained 1,868,388 events from the TAMA300 
data and 1,292,630 events from the LISM data. For these events, we perform the time, mass and 
amplitude selections discussed in the previous section. In Fig. [51 we show a scatter plot of the events 
after coincidence selections in terms of p TAMA and Plism- A significant number of events are removed by 
imposing coincidence conditions. Only 0.04% of the TAMA300 events remain. In Table ITVl we show the 
number of events which survived after the selections. 

We reduce the fake events by introducing the renormalization p by \ 2 m addition to the coincidence 
conditions. In Fig.EH we show a scatter plot of these events in terms of the value of Ptama/ y Xtama an d 

Plism/ Valism- 

In order to obtain statistical significance from the above results, the number of coincident events should 
be compared with the number of accidental coincidences produced purely by noise events. If events occur 
completely randomly, and its event rate in each detector is stationary, the average number of accidental 
coincidences after the time selection is given by 

window 

Hpr = N TAMA A^ LISM ^ j (18) 

-L obs 

where A TAMA and A LISM are the number of events in each detector, T b s is the total observation time, and 
^window - 1S ^ e averaged value of the time selection window. The averaged value of the time selection 
window is evaluated as Ai™ ndow = 1.29 msec. We thus obtain n pl — 6.3 x 10 3 , which is slightly larger 
than the observed number of coincidence, 4706, after the time selection. One reason for this diffrence is 
that the event trigger rate is not stationary over the whole period of this observation. 

In order to obtain a more reliable value for the rate of accidental coincidence, we use the time shift 
procedure. Namely, we shift all events of one detector by a time 8t artificially (which is called the time 
delay), and perform coincidence searches to determine the number of accidental events n c (6t) for various 
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FIG. 10: (ptama/V Xtama ,Plism/\/ Xlism) scatter plots. The crosses (+) are the events survived after the time 
selection, and the circled crosses (©) are the events survived after the time, mass and amplitude selections. 



Results of independent matched filtering searches 





TAMA300 


LISM 


Number of events 


1,868,388 


1,292,630 


Results of coincidence an 


alysis 








ft-acc i 0"acc 


after time selection 
after time and mass selection 
after time, mass and amplitude selection 


4706 
804 
761 


(4.2 ± 0.5) x 10 3 
(7.1 ±0.8) x 10 2 
(6.7 ±0.8) x 10 2 


Threshold 


iVobs 




Ptama/\/ Xtama > 8.3 and Plism/ \J Xlism > 8.1 





0.063 



TABLE IV: Results of coincidence analysis. n b s is the number of coincidence events. n acc and <r acc are the esti- 
mated number of accidental coincidence and its variance, respectively. Note that the mean number of accidentals 
and their variance after the time selection procedure affect those after the time and mass selection procedure, 
and the latter affect those after the time, mass and amplitude selection procedure. Thus, because the observed 
number of coincidence events is consistent with the expected number of accidental coincidence after the time 
selection procedure, it is not unnatural to find a good agreement between the observed value and the expectation 
value in each of the subsequent selection procedures. 



values of St 0] . With m different values of time delay, we calculate the expected number of coincident 
events and its standard deviation as 



711 < J 



(19) 
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FIG. 11: From top to bottom, the time delay histogram after time selection, after time and mass selection, and 
after time, mass and amplitude selection, respectively, are plotted. 



\ 



i=l 



(n c (SU) - n acc ^j /(m-1). 



(20) 



Since there is no real coincidence if \5t\ S> Aidis, the distribution of the number of coincidences with 
time delay can be considered as an estimation of the distribution of accidental coincidences. The number 
of coincident events, n c (0), is compared to the estimated distribution. 

Fig. II II shows the time delay histograms with m — 400. The 400 time delays are chosen from —12000 
sec to 12000 sec in increments of 60 seconds. The distribution of accidentals is shown in Fig. In 
Table Hvl we also list the expectation values of the number of accidental coincidence and the standard 
deviation after each selection procedure. As can be seen from this, the number of coincident events after 
each selection procedures is consistent with the expected number of accidental coincidences within the 
statistical fluctuations. Thus, we conclude that no statistically significant signals of real coincident events 
are observed in our search. 



VI. UPPER LIMIT TO THE EVENT RATE FROM COINCIDENCE ANALYSIS 

In this section, we present a method to evaluate the upper limit to the event rate based on the above 
result of the coincidence analysis. 

The upper limit to the event rate is given by Eq. (|10f> as in the case of the single-detector searches. The 
upper limit N to the average number of real events can be determined by Eq. (|12H . using the observed 
number of events N Q ^ S with pj \px^ greater than the threshold, the estimated number of fake events N^ g 
with pi \px* greater than the threshold, and the confidence level. We set different thresholds to the value 
of Ptama/ \J Xtama an( i Plism/ \J Xlism respectively. An advantage of this is that, because of its simplicity, 
it can be readily applied to the cases when more than two detectors with different arm directions are 
involved. 

We determine a background distribution f(y\, y 2 ) of the number of coincident events from the data for 
yi > 5.5 or y 2 > 5.5 in Fig. [TOI where yi = Ptama/\/Xtama an( i 2/2 = Plism/VxLsm- We evaluate the 
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FIG. 12: The distribution of the number of realizations with 400 different time delays with respect to the number 
of coincidences derived from Fig. 1111 From top to bottom, the distribution after time selection, after time and 
mass selection, and and after time, mass and amplitude selection are plotted. 



expected number of fake events which p/yx 2 is greater than the thresholds y\ = y T or y 2 = y L by 



/>oo />oo poo poo />oo />oo 

N bg = J dyi / dy 2 f(y 1 ,y 2 ) + / dyi J dy 2 f(y x ,y 2 )- / dyi / dy 2 f(yi,y 2 ). (21) 
J vt Jo Jo J y-L J w J v-l 



As the false alarm rate, we adopt 0.00023 [1/hour] (= 2.0 [1/yr]) which corresponds to the number 
of expected fake events Ab g = 0.063. We choose the thresholds y± = y T = 8.3 for TAMA300 and 
2/2 = 2/l = 8.1 for LISM. The observed number of events with y\ or y 2 greater than the threshold is 
A'obs = 0. Therefore we obtain the upper limit to the average number of real events with y\ or y 2 greater 
than the threshold as N = 2.30 (CL. = 90%) from Eq. |Q . 

The detection probability e is derived by the method explained in Section fV Bl and is shown in Fig. 1131 
With the thresholds chosen above, we obtain e = 0.182. Using the upper limit to the average number of 
real events N with y\ or y 2 greater than the threshold, the detection probability e and the length of data 
T = 275 [hours], we obtain an upper limit to the event rate within lkpc to be N/(Te) — 0.046 [1/hour] 
(CL.= 90%). 

Unfortunately, this value is not improved from the value obtained by the analysis of the TAMA300 
data. The dominant effect that causes the difference in the upper limit for a single-detector analysis and 
the coincidence analysis comes from the difference in the detection efficiency. The detection efficiency of 
the coincidence analysis in our case is determined by that of LISM, since LISM has the lower sensitivity. 
The efficiency of LISM is improved in the case of the coincidence analysis, since the threshold is lowered. 
However, this does not compensate the difference in the detection efficiency between TAMA300 and 
LISM. The efficiency of TAMA300 is already nearly 100 % in lkpc without performing the coincidence 
analysis. Thus, by taking the coincidence with the detector which has much lower sensitivity, the detection 
efficiency of the coincidence analysis becomes lower than the case of TAMA300 alone. As a result, the 
upper limit to the event rate we obtained by the coincidence analysis is less stringent than the one 
obtained by the analysis of the TAMA300 data. 
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FIG. 13: Detection efficiency in the coincidence analysis of sources within lkpc as a function of TAMA300 and 
LISM thresholds for p/^/x*. 



VII. SUMMARY AND DISCUSSION 



In this paper, we performed a coincidence analysis using the data of TAMA300 and LISM taken during 
DT6 observation in 2001. 

We analyzed the data from each detector by matched filtering and obtained event lists. Each event 
in the lists was characterized by the time of coalescence, masses of the two stars, and the amplitude of 
events. If any of the events are true gravitational wave events, they should have the consistent values 
of these parameters in the both event lists. We proposed a method to set coincidence conditions for 
the source parameters such like the time of coalescence, chirp mass, reduced mass, and the amplitude of 
events. We took account of the time delay due to the distance between the two detectors, the finite mesh 
size of the mass parameter space, the difference in the signal amplitudes due to the different sensitivities 
and antenna patterns of the detectors, and errors in the estimated parameters due to the instrumental 
noise. Our Monte Carlo studies showed that we would not lose events significantly by imposing the 
coincidence conditions. 

By applying the above method of the coincidence analysis to the event lists of TAMA300 and LISM, 
we can reduce the number of fake events by a factor 10~ 4 compared with the number of fake events 
before the coincidence analysis. In order to estimate the number of accidental coincidences produced by 
noise, we used the time shift procedure. We found that the number of events survived after imposing 
the coincidence conditions is consistent with the expected number of accidental coincidences within 
the statistical fluctuations. Thus we found no evidence of gravitational wave signals. As discussed in 
AppendixJCJ the sidereal time distribution of the survived events were also consistent with the distribution 
of accidentals. 

Finally, we proposed a simple method to set an upper limit to the event rate and applied it to the 
above results of the coincidence analysis. We obtained an upper limit to the Galactic event rate within 
lkpc from the Earth to be 0.046 [1/hour] (90% CL). In our case, since LISM has a much lower sensitivity 
than TAMA300, we were unable to obtain a more stringent upper limit to the event rate than the one 
obtained by the single-detector analysis of TAMA300. This is because the detection efficiency in the 
coincidence analysis is determined by the detector with a lower sensitivity. 

However, if we have two detectors that have comparable sensitivities, it is possible to obtain an improved 
upper limit compared to a single-detector analysis. As an example, let us imagine the case when the 
sensitivity of LISM is the same as that of TAMA300. The result of Galactic event simulations suggests 
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that the detection efficiency in the case of a single-detector analysis is 0.35, while it improves to 0.48 in 
the case of a coincidence analysis. These values are translated to upper limits on the Galactic event rate 
of 0.026 [1/hour] (90% CL) for the single-detector case and 0.019 [1/hour] (90% CL) for the two-detector 
case. 

The method of a coincidence analysis and the method to set an upper limit to the event rate proposed 
here can be readily applied to the case when there are more than two detectors with arbitrary arm direc- 
tions. Hence these methods will be useful for data analysis for a network of interferometeric gravitational 
wave detectors in the near future. 
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APPENDIX A: METHOD TO DISTINGUISH BETWEEN REAL EVENTS AND 

NON-GAUSSIAN NOISE 

The real data from TAMA300 and LISM contain non-stationary and non-Gaussian noise. One way to 
remove the influence of such noise is a veto analysis by using the data of various channels which monitor 
the status of the interferometers and their environments. Such an analysis has been performed using the 
data of TAMA300 24] . However, more efforts will be needed to establish an efficient and faithful veto 
method. 

It was shown that about 20 % of the data from TAMA300 DT6 contains non-Gaussian noise signifi- 
cantly 20]. Even if we remove this portion of the data with large non-Gaussian noise, the rest of data 
may still contain some non-Gaussian noise. It is thus necessary to introduce a method by which we can 
discriminate the non-Gaussian noise from real gravitational wave signals using the properties of inspiral 
signals. As one of such methods, the % 2 method was introduced in 0]. 

In this method, we examine whether the time-frequency behavior of the data is consistent with the 
expected signal. We divide each template into n mutually independent pieces in the frequency domain, 
chosen so that the expected contribution to p from each frequency band is equal: 

kcMf) = mUif) + *£.,(/) + ■ • • + %1,C/)- (Ai) 

We introduce 

*£.) = M2.)>' *(t) = ~( s '^))- (A2) 

Then, \ 2 is defined by 




with 

«io = (^.*S§) = (C fc (S) = i (A4) 

Provided that the noise is Gaussian, this quantity must satisfy the x 2 -statistics with 2n — 2 degrees of 
freedom and is independent of p — sl^f^ + z f s ) ■ F° r convenience, we use a reduced chi-square defined by 

X 2 /(2n — 2). In this paper, we choose n = 16. 

In the case of TAMA300, it was found that there was a strong tendency that noise events with large 
X 2 have large values of p. Since the value of \ 2 will be independent of the amplitude of inspiral signals 
when the parameters such as t c , M and rj of the signal are equal to those of a template [2jJ, one may 
expect that we can discriminate real signals from noise events by rejecting events with large x 2 > and this 
method was used in the TAMA300 DT2 analysis 0. 

However, in reality, since we perform analysis on a discrete t c and a discrete mass parameter space, 
the parameters of a signal do not coincide with those of a template in general. We have found in the 
analysis of the TAMA300 DT4 data in 2000 that this difference produces a large value of \ 2 when the 
SNR of an event is very large even if the event is real Thus, if we apply a threshold to the value 
of x 2 to reject noise events, we may lose real events with large SNR. This is a serious problem since an 
event with a large SNR has a high statistical significance of it to be real. This lead us to introduce a 
different rejection criterion when we performed an inspiraling wave search with the TAMA300 DT4 data 
|22j . namely, a threshold on the value of p/ y/x 2 - By Galactic event simulations, we found that this new 
criterion can give a better detection efficiency of the Galactic events without losing strong amplitude 
events. 

Here we examine whether the pj \J~x* selection is useful also in the case of the TAMA300 DT6 data. 
For comparison, the detection efficiency for a simple \ 2 threshold is shown in Fig. 1141 
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FIG. 14: Detection efficiency for a x 2 threshold. In this figure, the threshold is set to \ 2 = 1-5. 

For the x 2 < 1-5 threshold, using 287.6 hours of the data, the false alarm rate 2.0 [1/yr] determines 
the SNR threshold to be p = 12.5. This gives the detection efficiency of 0.213. On the other hand, 
as discussed in Section llYl the detection efficiency in the case of the pf \/~x 2 threshold is 0.263 for the 
same false alarm rate, 2.0 [1/yr]. We thus find that we have a better efficiency for the p/-\/a threshold, 
although the gain of efficiency is not very large. However, the important point is that we have much 
larger detection efficiency for signals with large SNR. 



APPENDIX B: DIFFERENT CHOICE OF At 

In this appendix, we consider the case of a different choice of the length of duration At to find local 
maximum of matched filtering output (see IIIIB|l . to see if our conclusion is affected by a different choice 
of At. 

Here we adopt At = 3.28sec. In this case, the total number of events is found to be 158,437 for 
TAMA300 and 142,465 for LISM. The numbers of events survived after each step of the coincidence 
selections are given in Table Ivl The corresponding estimated numbers of accidentals are also shown. The 
scatter plots of these selected events are shown in Figs.ll5landll6l We see that the number of coincident 
events is consistent with the number of accidentals within the standard deviation, in agreement with our 
conclusion given in the main text of this paper. 



APPENDIX C: SIDEREAL TIME DISTRIBUTION 

In this appendix, we examine the sidereal time distribution of the events. In Fig. I17l (ah we plot the 
number of coincident events as a function of the local sidereal hour at the location of TAMA300. The 
estimated number of accidental coincidences are also plotted, which are obtained by the same time shift 
method used in Section IV CI but for data within each bin of the sidereal hour. If the gravitational wave 
sources are sharply concentrated in the Galactic disk, we would detect more events when the zenith 
direction of the detector coincides with the direction to the Galactic plane than the rest of time. The 
zenith direction faces to the Galactic disk at around 6:00 and 18:00 in the sidereal hour. Since LISM is 
only sensitive to sources within a few kpc, we may not be able to see any significant excess of the events in 
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FIG. 15: (ptama,Plism) scatter plots in the case At = 3.28sec. The crosses (+) are the events survived after the 
time selection, and the circled crosses (©) are the events survived after the time, mass and amplitude selections. 
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FIG. 16: (p tama/v Xtama jPlism/y Xlism) scatter plots in the case At — 3.28sec, The crosses (+) are the events 
survived after the time selection, and the circled crosses (©) are the events survived after the time, mass and 
amplitude selections. 



the Galactic disk within this distance unless the concentration of the sources to the Galactic disk is very 
strong. Even in this case, it is useful to investigate the sidereal time distribution to look for signatures 
of real events. 

We find that the distribution of coincident events is consistent with accidentals, although there are a 
few hours in which the agreement is not very good. Thus, we conclude that the result of the sidereal 
hour distribution is consistent with the number of accidentals, and there is no signature of gravitational 
wave event. 

In Fig. 1171 (b). we also plot the number of coincident events as a function of the Japanese Standard 
Time ( JST) . Since the deviation of the local sidereal time from JST is not very large during the period 
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Results of independent matched filtering searches 

TAMA300 LISM 
Number of events 158,437 142,465 



Results of coincidence analysis 



^obs ^acc i (7 a 



after time selection 70 75.0 ± 8.6 

after time and mass selection 18 18.8 ± 4.1 

after time, mass and amplitude selection 17 17.9 ± 3.8 



TABLE V: Results of coincidence analysis in the case At = 3.28sec. n b s is the number of coincidence events. 
n acc , Oacc are the estimated number of accidental coincidence and its variance. 




4 8 12 16 20 
Local sidereal time [hour] 
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12 16 
JST [hour] 



FIG. 17: (a) The event distribution as a function of the local sidereal time. The solid line represents the number 
of coincident events per one sidereal hour. The dot-dashed line represents the estimated number of accidentals, 
(b) The number of coincident events as a function of the Japanese Standard Time. 



of observation, this figure is very similar to Fig.EI( a )- The reason that there are many coincident events 
during 20:00 to 22:00 JST is due to a large number of events recorded by LISM during that period. 
During the DT6 observation, there were some activities in the Kamioka mine from 20:00 to 22:00 JST, 
and trucks went through the tunnel of the mine during that period. We suspect this caused fake events 
in LISM. 
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APPENDIX D: PARAMETER ESTIMATION ERRORS INDUCED BY DETECTOR NOISE 



In this appendix, we briefly review the theory of the parameter estimation error developed in [23j . This 
is used in determining the parameter windows for the coincidence analysis in this paper. 

In the matched filtering, for a given incident gravitational wave, different realizations of the noise will 
give rise to somewhat different best-fit parameters. For a large SNR, the best-fit parameters will have 
Gaussian distributions centered on the correct values. Specifically, let 9 l be the correct values of the 
parameters, and let 9 % + A9 l be the best-fit parameters in the presence of a realization of noise. Then 
for large SNR, the parameter estimation errors A9 Z have the Gaussian probability distribution 

p(A0 i )=JVe-* r « AO<A * i . (Dl) 
where F^j is the called Fisher Information matrix defined by 

/ dh dh\ , . 

V ^\m-W& (D2) 



and M — ^/det(r/27r) is the normalization factor. It follows that the root-mean-square errors in 9 l is 
given by 



a. t = v/((A^) 2 ) = VE^, (D3) 
where S = and the correlation coefficient between parameters 9 l and 9^ is given by 

„ < A'/'A'/-' > £y 



(D4) 



By definition, each c y lies in the range (—1,1). 

As given in Section IlII Al an inspiraling signal in the frequency domain is given by 

h(f)=Af- J/6 e l ^ f) . (D5) 

Here we consider the phase tp(f) up only to the second post-Newtonian order but including the effect of 
the spins of stars. Note that this is slightly different from the template formula (JIJ used in our analysis. 
The phase ip(f) is given by 



^(/) = 2-Kft c -4> c - J + I ^( 7 r>l/)- 5 / 3 [l + ^(^ + ^ ?7 )(7rM/) 2 / 3 -4(47r-/3)(7rM/) 



20 /743 11 



'3058673 5429 617 
,1016064 + 1008 77 + 144 

In the above, (3 is the spin-orbit parameter given by 



/30586Y3 542y 617 , \. „,„, 4 /ol /T ^„, 

10 h r? H 77 2 - cr ttM/ 4/3 . (D6 

V 1016064 1008' 144' K J1 J K 1 



1 2 



12 

i=l 



and ^ i = Si/mf, and Si is the spin angular momentum of each star, and L is the unit vector along the 
orbital angular momentum vector. The spin-spin parameter a is given by 



We define 



<7=^(-247x 1 -X 2 + 721L. X iL-x 2 ). ( D8 ) 



f-7/3 

f=**l L—tf. (D9) 
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Wc also define the frequency moments f a of the noise spectrum density: 

/r/a = df [f^SMT 1 (D10) 

Jo 

U = fyl df {rSMT 1 - (Dll) 

In order to evaluate the Fisher matrix, we calculate the derivatives of h(f) with respect to the seven 
parameters 

= QnA,f o t e ,<l>c,hiM,biTi,P,<T), (D12) 

where /o is a fiducial frequency which is taken to be the frequency at which S n (f) becomes minimum. 
We obtain 

al{1) - Sen, 



d\nA 

i~Kf), 



dh(f) „ / / 
7777^ = 2m 

dh(f) 



dh(f) = _^ r? -4/5 (7rM/) - 1 /3^ (/) _ (D13) 



9cr 64 
Here we have defined 



and 



4 /743 11 
4 ~ 3V336 + T ?/ 

B 4 = |(47r-/3), 
5 

/ 3058673 5429 617 2 \ 

C 4 = 2 1 77 H 77 -tr , (D14) 

V 1016064 1008' 144' 7' v ; 



743 33 

^ = 168-T 77 ' 
27 

5 

/3058673 5429 617 2 
Cs = 1 ( 1016064 + 4032 ^ + ~96~ ~ 



Finally, the components of T can be obtained by evaluating Eq. i|D2l) . They can be expressed in terms 
of the parameters 9, the signal-to- noise ratio p, and the frequency moments fg. The components of T^- 
are given by 

TinAj = SinAjP 2 , {j = tn.A,f t c , <p c ,]nM,bir),o;0), (D16) 
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r tc t c = (2^) 2 |/" 1/3 p\ (D17) 
Jo 

T tc c = -27r^/ 4/3 p 2 , (D18) 
Jo 

r tc i„,M = -^(^)" 6/3 (/3 + ^(ttM) 2 / 3 - B 4 (7rM)/ 2 + C 4 (ttM) 4 / 3 / 5/3 ) p 2 , (D19) 

rt e in, - -^(^)- 5/3 (A5(7rM) 2 / 3 -B5(^M)/ 2 + C5(^M) 4 / 3 / 5/3 ) p 2 , (D20) 

I\ ln/3 = -^ v -^(nM)-^f 2 p\ (D21) 

IV = -^f ^ 4/5 (--M)- 1/3 /v3 ^ (D22) 
-J A/o 

r 0ctc = r tc0c , (D23) 

1^0. = p 2 , (D24) 
r 0c inA< = ^(7rA1)- 5 / 3 (/4 + A4(7rM) 2 / 3 / 10/ 3-i? 4 ( 7 rM)/3 + C4(^M) 4 / 3 / 8/ 3) p 2 , (D25) 

r^in, - ^(^)- 5 / 3 (A 5 (^M) 2 / 3 A 0/3 -B 5 (7rM)/ 3 + C 5 (7rM) 4 / 3 / 8/3 ) p 2 , (D26) 

T, c , = -l„-3/6(7T>W)- 2 /3/3 p\ (D27) 

T, cCT - g^ 4/5 (--M)- 1/3 / 8 /3 P 2 , (D28) 
TinTWtc = r tclnj vi, (D29) 

TlnA^^c = T^inX, (D30) 
TlnMlnM = 16384( f^ ) 10/3 (^(^ M ) 4/3 /l7/3/l3/3 + 2A 4 (ttM) 2 / 3 {/ 15/3 - B 4 nMf 4 

+ C 4 (^M) 4 / 3 / 11/3 + irM(BlirMf 11/3 - 2B 4 (f 14/3 

+ C 4 (nM) 4 ^f 10/3 ))}) p\ (D31) 

rin - Mln " = 12288^8/3^10/3 (5M 2 / 3 (A 5 {/ 15/3 + A 4 (nM)^ f 13/3 - B 4 nMf 4 
+ C 4 (7rM) 4 / 3 / 11/3 + (nM)^ 3 (C 6 {nM)^ 3 {f 13/3 + A 4 (nM) 2 / 3 f u/3 

- B 4 nMf 10/3 + C 4 (irM)^ 3 f 3 ) - B 5 (/ii4/3 + A 4 (^M) 2 / 3 / 4 

- B 4 7rMf u/3 + C 4 (7rM) 4 / 2 / 10/3 )))}) p 2 , (D32) 

r lnM p = -^V- 3/5 (*M)- 7 ' 3 (f 1A/3 + A 4 {^Mf/ 3 j 4 - B 4 (7rM)/ 11/3 

+ C 4 (^M) 4 / 3 / 10/3 ) p 2 , (D33) 

7^ / - 

rm^. = 8192^ 4/5 (^) _2 (/i3/3 + ^M) 2 / 3 / 11/3 - B 4 {-kM)J w/3 

+ C 4 (7rM) 4 / 3 / 3 ) p 2 , (D34) 
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■Tin rjt c — ~^t c In 77 1 

(D35) 

rinr;0 c = ^4> a \nrn (D36) 
rinr)lnX = "F\nM\nni (D37) 

Tin, i„, = ^xe^iQ/a (^ 4/3 {^/"i3/3 " lA^M) 1 '* f A + S 5 2 (ttM) 2 / 3 Al/3 

+ 2A 5 B 5 (irM) 2 / 3 f 11/3 - 2B 5 C 5 TrMf 10/3 + C 5 2 (^A/) 4 / 3 / 3 }) p 2 , (D38) 

IW = -^^- 3/5 (7rA4)- 7 / 3 (A( 7 rM) 2 / 3 /4-B 5 (7rM)/ n/ 3 

+ C 5 (7rAf) 4 / 3 /io/3) P 2 > (D39) 
IW = ^v- i/5 ^Mr 2 (A 5 (7TM) 2 / 3 f 11/3 -B 5 (7rM)f w/3 

+ C b {TiM) A / z f 3 ) p 2 , (D40) 



I>t c = r te/J , (D41) 

r Wc = r 0c/3 , (D42) 

T/31nA4 = Fin MP > (D43) 
T/Jlnr) = ri n??/ 3, (D44) 

^ = l^i^" 675 ^)^ 3711 / 3 ^' (D45) 

IV = -^^ 7/5 (vrM)- 1 / 10 / 3 P\ (D46) 

r fftc = r tcCT , (D47) 
r CT0c = r 0cCT , (D48) 

r<jln.M = Tin^vicr, (D49) 
rcrln?7 = ^Intyai (D50) 

r ff/3 = IV, (D5i) 
r ™ = S^ 8/5( ^ r2/3/3 p2 - (D52) 

It is ensured by these formulas that the eigenvalues of the Fisher matrix are always positive definite. 

The variance-covariance matrix can now be obtained from S = T -1 , and the root-mean square 
errors and the correlation coefficients are computed from Eqs. I|D3|I and l|D4(l . 

For example, using a typical noise spectrum density of TAMA300, the root-mean square errors of 
the parameters in the case p = 10 and (3 = a = are evaluated to be Ay4. TAMA /^4 TAMA = 0.10, 
AtJ AMA = 0.65msec, A<^ AMA = 6.88radians, A7W TAMA /7W TAMA = 1.43x lO" 2 , and Arj TAMA /r] TAMA = 
2.47 x 10" 1 . 
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